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Ch This work is a contribution to the study of universality in out-of-equilibrium 

I ' lattice models undergoing a second-order phase transition at equilibrium. The 

"""O ■ experimental protocol that we have chosen is the following: the system is pre- 

£h ' pared in its high-temperature phase and then quenched at the critical temper- 

ature T c . We investigated by mean of Monte Carlo simulations two quantities 
that are believed to take universal values: the exponent A/z obtained from 
the decay of autocorrelation functions and the asymptotic value Xoo of the 
£NJ , fluctuation-dissipation ratio X(t,s). This protocol was applied to the Ising 

^ ■ model, the 3-state clock model and the 4-state Potts model on square, tri- 

r**"* ' angular and honeycomb lattices and to the Ashkin- Teller model at the point 

belonging at equilibrium to the 3-state Potts model universality class and to a 
■ multispin Ising model and the Baxter- Wu model both belonging to the 4-state 

Potts model universality class at equilibrium. 

1 Introduction 



Abstract 



Universality is an extremely fruitful concept in statistical physics and has been 
widely studied in the context of systems undergoing a second order phase transition. 
At thermodynamic equilibrium, the length scale £ of spatial correlation functions 
of the local order parameter diverges as a power-law with a critical exponent v as 
the temperature approaches the critical temperature. Other observables (energy, 
magnetisation, . . .) display such a power-law dependence as well but with different 
critical exponents. It turns out that many microscopic details of the Hamiltonian 
do not change the value of these critical exponents. The big success of the renor- 
malisation group has been to explain that a few number of these critical exponents 
are independent and that different models have the same set of critical exponents if 
they differ only by irrelevant operators. Consequently, the usual models of statisti- 
cal physics can be classified into universality classes according to the value of their 
critical exponents. The space and order parameter dimensions, the Hamiltonian 
symmetries, the presence of long-range interactions, of randomness . . . determine 
the universality class. Not only critical exponents but also ratios of critical ampli- 
tudes turn out to be universal. 

The question of universality in out-of-equilibrium processes has been addressed 
in the context of dynamical transitions undergone by reaction-diffusion systems [1]. 
A set of exponents can be defined, for instance using the algebraic decay with time 
of the density of active sites. 
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In the following, we will restrict ourselves to another important set of out- 
of-equilibrium processes provided by systems that undergo a second-order phase 
transition at equilibrium and that are quenched from their high-temperature phase 
to their critical point or their low-temperature phase. Because of the competition 
of domains in different equilibrium states, such systems never reach equilibrium in 
the thermodynamic limit [2]. The hypothesis was made that systems belonging to 
the same universality class at equilibrium share universal quantities if their time- 
evolution is governed by dynamics satisfying the same conservation laws and if order 
parameter and conserved quantities are related in the same way [3] . In the following, 
we will consider systems belonging to the same universality class at equilibrium and 
whose dynamics has no conserved quantity (model A). In this context new exponents 
have been defined: the dynamical exponent z related to the growth of the domain 
length scale L{t) with time by [3] 

L(t)~t 1/Z (1) 

and the autocorrelation exponent A related to the decay of two-time autocorrelation 
functions of the local order parameter [4, 5] by 

C{t,s) ~ t- x ' z . (2) 

These two new exponents z and A take different values whether the system is 
quenched at its critical temperature T c or below. They are believed to be universal. 
Numerical calculations for the Ising model on square, triangular and honeycomb 
lattices [6] or for three different models belonging to the Ising equilibrium univer- 
sality class [7] support this conjecture for the exponent z aXT ~ T c . Simulations for 
other models, for instance the 3-state Potts model [8], give estimate for z sufficiently 
close to give support to the conjecture of a dynamic universality. 

The decay of persistence, i.e. the probability that the total order parameter 
has not changed sign at time t, defines a third new non-trivial exponent 9' . In the 
rather unusual case where the dynamics of the total order parameter is Markovian, 
it is related to the two previously-defined ones by z6' = A — d + 1 — r]/2 [9]. 

Renormalisation-group study of the 0(n) model in d = 4 — e dimension shows 
that if the system is quenched at the critical temperature from an initial state with 
a small non-vanishing magnetisation, the latter first grows as t e before decaying 
asymptotically as ^ fj/lJ ~ L(t)- & l v ~ t-< 3 / l ' z [10]. The initial-slip critical exponent 
9 is related to the autocorrelation exponent by 

\ = d-z6 (3) 

where d is the dimension of the space. This behaviour of the magnetisation has 
been exploited in the so-called short-time dynamics Monte Carlo method [11]. The 
question of universality has been addressed by several authors in this context. The 
Ashkin- Teller model has been studied by short-time dynamics for several points 
of its exactly-known critical line [12]. Unfortunately, the point belonging to the 
3-state Potts model universality class, i.e. y — 3/4 where v — jE^, has not been 
considered. Noting that the estimates of 9 given by the authors vary roughly linearly 
with the parameter y, one can estimate the exponent 9 at the point y = 3/4 to be 
9 ~ 0.111. This result is not compatible with that obtained for the 3-state Potts 
model [8]: 9 ~ 0.0815(27). Moreover, numerical estimates of this exponent 9 for 
the Baxter- Wu [13] and a multispin Ising model [14] both belonging to the 4-state 
Potts model equilibrium universality class are incompatible with estimates for the 
latter while the dynamical exponent z seems to be the same [13]. 

Apart from exponents, ratios of amplitudes have also been conjectured to be 
universal, as for example the fluctuation-dissipation ratio X(t, s) in the asymptotic 
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limit t, s — ► +00. The response function R(t, s) to an infinitesimal field h(s) coupled 
to the local order parameter is expected to decay algebraically with time with the 
same exponent X/z as autocorrelation functions provided the initial state does not 
display spatial correlations. Note that autocorrelation and response functions are 
related at equilibrium by the fluctuation-dissipation theorem (FDT): 

k B TR(t,s) = -^- s C(t,s). (4) 

On the basis of a mean-field study of a spin-glass-like model, it was conjectured 
that the fluctuation-dissipation theorem may be generalised by adding to Eq. 4 a 
multiplicative factor depending on time only through autocorrelation functions [15]: 

k B TR(t,8)~X[C(t,8)]^C(t,s), (t~«»l). (5) 

However, numerical simulations for the Ising model [16] and renormalisation-group 
calculations for the 0{n) model [17] suggest that the fluctuation-dissipation ratio 
is not a function of C(t,s) only but of t/s. Scaling arguments constrain the au- 
tocorrelation and response function to the following asymptotic behaviour [18] at 

C{t,s)~s-%f c (^), k B T c R(t,s)~s- 1 -%f R (^) (i~s»l) (6) 

where both fc{x) and fn{x) are scaling functions whose asymptotic behaviour is 
given by 

fc/ R {x) ~ A c/R x- x/z , x -» +00. (7) 

Inserting Eq. 6 and 7 into Eq. 5 leads to the asymptotic value of the fluctuation- 
dissipation ratio: 

(8) 



A R 

Xqo = lim X(t, s) = — — 
M^+°o A c 



X 2/3 

z vz 



which turns out to depend only on exponents which are believed to be univer- 
sal and on the ratio of autocorrelation and response amplitudes. The latter has 
been conjectured to be universal [18]. This conjecture applies thus to as well. 
Renormalisation-group calculations [17] of the O(l) model gives indeed an estimate 
of Xoo compatible with numerical values for the Ising model [19, 16, 21]. Moreover, 
a recent numerical calculation of the integrated response function to an exchange 
coupling perturbation gave an estimate of X^ in full agreement with theses val- 
ues [19]. The Ising model may be peculiar since a one-loop renormalisation-group 
calculation of the O(l) model gives the same value of X^ whether the order pa- 
rameter is coupled to a conserved quantity (model C) or not (model A) [20]. Let 
us mention that calculations for the Ising-Glauber chain indicates that X^ takes 
a non-vanishing value only with a well defined protocol : a quench at the critical 
temperature from an initial state with an infinite number of domain walls, i.e. from 
an initial disordered state [22]. 

The plan of this paper is the following: in the first section, the expressions of the 
Hamiltonians of the different models we studied are given. The Glauber dynamics 
is then defined and we review the method we used to calculate the fluctuation- 
dissipation ratio without resorting to the Cugliandolo conjecture (Eq 5). The sec- 
ond section is devoted to the characterisation of the three universality classes under 
consideration. The question of the influence of the lattice on the quantities that are 
supposed to be universal is addressed in the third section. The fourth section is de- 
voted to the comparison of these universal quantities for different models belonging 
to the same universality class at equilibrium. 
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2 Dynamics and Models 



2.1 Definition of the models 

We have considered two-dimensional lattice models belonging at equilibrium to 
three different universality classes. The first of them is the Ising model [23] defined 
by the Hamiltonian: 

Rising = - J y^g'i^j ~ y]hjVi> Cj = ±1 (9) 

where the first sum extends to nearest neighbours on the lattice. A local magnetic 
field coupled to the local order parameter m t = Ui has been added to allow for the 
definition of a response function. For vanishing magnetic field, the model undergoes 
a second-order phase transition associated to the breaking of a Z(2)-symmetry. 
Duality relation leads to the critical point on the square lattice (3 C J — \ ln(l + \/2) 
where j3 c — l/fcsT c [24]. Although triangular and honeycomb lattices are not self- 
dual but dual of each other, the exact determination of the critical temperature 
is nevertheless possible [25] and gives (3 C J ~ 0.2746 on a triangular lattice and 
(3 C J ~ 0.6585 on the honeycomb lattice. 

The second universality class we have considered is that of the 3-state Potts 
model [26] whose Hamiltonian is 

Wp otts = -J^bai^i ~ ^2 h i 5 <?ifl, (Ti = Q,...q-l (10) 

where the number of state q is set to q — 3. For vanishing magnetic field, the model 
undergoes a second-order phase transition associated to the breaking of a Z(3)- 
symmetry. Self-duality relation allows for the exact determination of the critical 
point: (3 C J = ln(l + ^fq) at zero-magnetic field on the square lattice. Critical 
points can also be obtained for the triangular and honeycomb lattices [25]. In 
order to achieve better numerical stability, we used an equivalent formulation of 
this model, known as the clock model, whose order parameter is vanishing at the 
critical temperature without resorting to any additional normalisation: 



Hciock = - 




= 0,...«-l. (11) 



The critical temperature is readily obtained to be two-third of that of the equivalent 
g-state Potts model. As a prototype of model belonging to this universality class, 
we have chosen the Ashkin- Teller model [27] defined by the Hamiltonian 

Hat = - J ^ crjo-j - J 1 y^ TjTj - K (Ti<JjTjTj — /^q-j, (7j,Tj = ±l. (12) 

(ij) (ij) (i,f) * 

Indeed, in the isotropic case J = J' the universality class changes from that of the 
Ising model to that of the 4-state Potts model along an exactly-known critical line. 
The point belonging to the 3-state Potts model is defined by the critical couplings: 
(3 C J = j3 c J' ~ 0.34763 and (3 C K ~ 0.14209 [28]. 

The third universality class we have considered is that of the 4-state Potts model. 
For vanishing magnetic field, the model undergoes a second-order phase transition 
associated to the breaking of a Z(4)-symmetry. Duality relations allows for the 
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determination of the critical point on the square, triangular and honeycomb lat- 
tices [25]. Apart from the 4-state Potts model, we have made calculations for two 
other models belonging to this universality class: the Baxter- Wu model [29] defined 
by the Hamiltonian 

H-BW = — J ^ GjQjCJk — hjUj, &i = ±1 (13) 

(i,j,k) i 

where the sum extend over all triangles of a triangular lattice and a multispin Ising 
model [30] defined on a square lattice by the Hamiltonian 

Hut = - J 0^ - K ^ <Tiaja k - ^ h^i, a l = ±1 (14) 

where the first sum extends over couples of neighbour spins along the horizon- 
tal direction and the second over triplets of spins aligned along the vertical direc- 
tion. For the Baxter- Wu model, the ground states correspond to spin configurations 
where triangle vertices are decorated with spins {+,+,+}, {+,—,—}, { — ,+,—} or 
{ — ,—,+}. In the case of the multispin Ising model, the ground-states correspond 
to rows decorated with the sequence of spins: {+, +, +}, {+, — , — }, { — , +, — } or 
{ — ,—,+}. Duality relations lead for the Baxter- Wu model to the critical point 
[3 C J ~ | ln(l + \/2) and for the multispin Ising model to the same critical line in 
the J — K plane than the anisotropic Ising model [30] from which we have chosen 
the point (3 C J = 1, C K ~ 0.1362. 



2.2 Glauber dynamics and observables 

The time evolution of all these models is governed by the same discrete-time Glauber 
dynamics. Let us denote by p({a},t) the probability of the spin configuration {a} 
at time t. The dynamics is defined by the usual master equation for Monte Carlo 
simulations: 

p(M, t + 1) = £ W({a'} {a}, t)p({a'}, t) (15) 
{<?'} 

where VF({er} — > {a'},t) is the transition rate per time step from the state {a} to 
the state {a'} at time t. The conditional probability p({er}, t|{<r'}, s) for the system 
to be in the state {a} at time t knowing that it was in the state {a 1 } at time s 
satisfies the master equation too. The condition of stationarity ^p oq ({o'}) = of 
the equilibrium Boltzmann probability peq({c}) is ensured by the so-called detailed 
balance condition 

P^iW})W{W} - {a}, t) = Pcq ({a})W({a} - {a'}, t). (16) 

This last unnecessary but sufficient condition is fulfilled by the single-spin flip dy- 
namics defined for the single-spin flip oi — > oi ^ cr^, where <7j is a randomly chosen 
spin, whose transition rates are 

w({* } ^ wu) = [n<w;] Sa>i ^;^ AE (it) 

where <7j is a trial state for the spin a t and AE = W(ai, . . . <7j, . . . cr N )—Ti(ai, . . . cr i} . . . a N ) 
is the energy difference when replacing <Ji by di . The Markov chain has to be aver- 
aged over all possible spin-flips to recover Glauber dynamics [31]. This is the usual 
way Monte Carlo simulations proceed. In the simulations to be presented in the 
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following, the sites on which spin flips are applied are randomly chosen. Let us 
mention that another definition of the transition rates may be given 

r-r-r i e-WiW}) 
W({a} -» {a%t) = ^ e - <w(ci ,... J< ,... w) (18) 

which corresponds to the Heat-Bath algorithm in Monte Carlo simulations. 

Autocorrelation functions of the local order parameter rrn have been defined as 
usual as 

Cij(t,s) = (mi(t)mj(a)) (19) 
= Yl m i ({a}M{a},t\W},8)m j {{ t /})p({a f },s) 

W>{"'} 

Since periodic boundary conditions are used, autocorrelation functions are expected 
to be invariant under space translations. In the following, we will consider aver- 
aged autocorrelation functions over the lattice: C(t,s) = ± J2i=i c u(t, s )- Thc 
response to an infinitesimal field hi coupled to the local order parameter has 
been computed using a recently proposed method for the Ising-Glauber model [16]. 
We shall derive this result for a more general spin lattice model. The coupling of 
thc infinitesimal field to the local order parameter is performed by modifying the 
transition rate Eq. 17 in the following way: 

This modification corresponds to the addition of thc Zccman Hamiltonian —hiirii 
to the Boltzmann weight of the equilibrium probability distribution in Eq. 16. Thc 
average order parameter at time t can be expanded as: 

( mj (t)) = X>i(M)p(M,*) (21) 
m 

= E rn 3 ({a})p({<j},t\{<j'}, S + l)p({a'}, S + l) 

M>{"'} 

= E ^ (M) ?(M, t\{</}, s + l)W({a"} {a'}, s)p({a"}, s). 
{»"} 

The magnetic field being branched only during the time step s, the transition rate 
VF({er"} — ► {c'}, s) is the only quantity depending on hi provided that the spin-flip 
occurring at time s affects the spin branched to the magnetic field. The derivative 
of the transition rate being 

= ,[ mi(M) _!-Mi±=M)^; 

where e^(s) is equal to 1 if the transition rate involves a spin-flip on the spin <jj at 
time s and otherwise, the response to the perturbation follows 

R ^ S) ^{^hW) = £ rn 3 ({a})p({a},t\{a'},.s + l) 

V ^ W h=0 {a}, {a 1 } 

dW 



ei (s)W({a"}^{a'},s) 



, , ( { a"}^{a'},s)ei(s)p({a"},s) 
K un t [S) j h=Q 

= f3(m j (t)[m i (s + l)-mj f (s)]ei(s)) (23) 
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where mj^(s) is the average local order parameter over <Ti(s) and the trial value bi 
that appeared in Eq. 22. The response is zero if the transition rate at time s involved 
any other spin than cr,. This procedure is easily implemented in Monte Carlo 
simulations. Since periodic boundary conditions are to be used, we will consider 
only averaged response functions over the lattice: R(t, s) = jj J2iLi Rii{t-> s )- Using 
the discrete-time master equation Eq 15 and the Boltzmann equilibrium probability 
distribution, the FDT can be shown to be kBTR(t,s) — C(t,s + 1) — C(t,s) at 
equilibrium for the above-define discrete-time process [16]. The ratio measuring the 
violation of the FDT out-of-equilibrium is readily obtained to be: 

k B TR(t,s) _ E t KfflM* + l)-mf( S )]e t ( S )) 
l ' SJ C{t,s + l)-C{t,s) £.< miW [ mi ( s + i)_ mi ( s )]) ' ^ > 

This method allows for the calculation of X(t, s) for all values of both s and t during 
the same Monte Carlo simulation without resorting to Cugliandolo conjecture Eq. 
5 and without applying any finite magnetic field that could induce non-linear re- 
sponse. This method has been applied to the calculation of the integrated response 
function [32] , the study of the XY- model [33] and of the Ising model with Kawasaki 
dynamics [34] . It turned out to require smaller values of s and t than measurements 
of the integrated response function to give an accurate estimation of X^ but on the 
other hand, an average over a larger number of histories is needed to give a stable 
estimate of X^ . 



3 Ising and Potts universality classes 

We first present the study of the Ising model, the 3-state clock model and the 4- 
state Potts model which belong to different universality classes at equilibrium. The 
system is first prepared at infinite-temperature and then quenched at its critical 
temperature T c at t = 0. Three lattices sizes were investigated: L = 128, 192 and 
256 in order to check the possibility of finite-size effects. Because of the particular 
structure of their ground-state, the multispin Ising model and the Baxter- Wu model 
were simulated on lattices of size L = 258, a multiple of 3. The data have been 
averaged over 10, 000 initial configurations for L = 128, 20, 000 for L = 192 and 
50, 000 for L = 256 and L — 258. These parameters are the same for all simulations 
presented in this paper. We measured the autocorrelation functions C(t, s) and the 
fluctuation-dissipation ratio X(t, s) for t < 1500 and s = 10, 20, 40, 80 and 160. For 
all observables, errors bars on the averages over the initial configurations have been 
estimated as the standard deviation. In principle, the response functions could be 
used to estimate the exponent X/z too but they turn out to be smaller and thus 
noisier than autocorrelation functions making the numerical estimates of exponents 
much less accurate. 



Data for the autocorrelation function C(t,s) have been grouped into bins of 
twenty points (forty for the 4-state Potts model whose fluctuations are much larger 
than other models). For each bin, an effective exponent X/z was measured by power- 
law interpolation over the points inside the bin. We took into account error bars 
by weighting each point with the inverse of its square error in the fit. The effective 
exponent can be considered local because values of t/s remain very close for points 
inside the same bin. The effective exponent is plotted on Figure 1 (on the left) versus 
the inverse mean position s jt of the bin. A fairly good collapse of curves correspond- 
ing to different values of s is observed for small values of t/s. The main difficulty 
in the determination of the asymptotic behaviour is that the effective exponent de- 
creases down to a plateau before slowly increasing as s/t is going to zero. Moreover, 
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Figure 1: Top: effective critical exponent X/z calculated by power-law interpolation 
of C(t,s) over the range t € [i m i n ; 1500] versus s/i m i n for the Ising model, the 3- 
state clock and the 4-state Potts model (from bottom to top). Bottom: fluctuation- 
dissipation ratio X(t, s) versus s/t for the same models (from bottom to top). The 
different symbols correspond to s = 10 (circle), s = 20 (square), s = 40 (diamond), 
s = 80 (triangle up) and s = 160 (triangle left). The lattice size is L — 256. 



autocorrelation functions are getting smaller and thus noisier when s/t is going to 
zero, adding to the difficulty Our final estimate for X/z is the intercept given by a 
linear fit over points in the range s/t £ [0; 0.05], corresponding to any value of s and 
with a weight corresponding to the inverse square error. The bold line on Figure 1 
corresponds to this fit. The drawback of this method is that the points with the 
smallest error bars, i.e. giving the largest contribution to the fit, correspond to the 
smallest values of s. As a consequence, problems may be caused by corrections to 
scaling depending on s only. However, as can be seen on Figure 1, the good collapse 
of curves corresponding to different values of s suggests that these corrections are 
weak. The data having been produced by Markov chains, the effective exponents, 
measured at different values oft or even s, are correlated. As a consequence, the lin- 
ear fit underestimates the true error. The assumption of an exponential decay of the 
autocorrelations G |^ = [j (t/s) — b — a|] ^(i'/s') — b — where a and b 
are the parameters given by a linear fit allows for the estimation of the correlation 
length I by summing the autocorrelations: I = J2t',s',t.s G (j' f) / Z)t,« G (§> f ) 



in the range s/t £ [0; 0.05]. The standard deviation on the intercept b as given by 
the linear fit is valid only in case of uncorrected data for which it decays as 1/ y/N 
where N is the number of degrees of freedom of the fit. To take into account the 
fact that only one point out of i is statistically independent, we multiplied the error 
on by \Tl. The correction is quite large since in the case of the Ising model for 
instance, the error is 3.10 -3 when neglecting correlations and 2.10 -2 when taking 
them into account. On Figure 1, our final estimate of the exponent X/z is given by 
the intercept of the bold line with the y-axis and its error bars are plotted along 
the y-axis. The values are collected in Table 1. They are compatible within error 
bars with the values found in the literature : X/z ~ 0.731(3) [35] or 0.732 [36] using 
z = 2.1667(5) [7] for the Ising model, 0.828(2) [8] for the 3-state Potts model and 
0.919 assuming z = 2.294 [13] for the 4-state Potts model. 

The scaling behaviour of C(t,s) and R(t,s), as given by Eq. 6, leads to a 
fluctuation-dissipation ratio X(t, s) depending only on t/s : 



k B TR(t,s) f R (f ) 

lC{t,s) ~ M /o (1) + i//, (i) 



JTCMH 'T V ~ m ,,V , 77V (25) 



Renormalisation-group calculations for the 0(n) model confirm that hypothesis [17]. 
Our numerical estimates of X(t, s) indeed collapse when plotted versus t/s. The 
asymptotic value = lim s ^ +00 lim t ^ +oc X(t, s) may thus be obtained as \\m. t / s ^ +OQ X{t, s). 
Note however that some exactly-solvable models display a cross-over when t becomes 
large so that the appropriate regime to take into account might not be s/t going 
to zero [19]. In our case, no significant dependence on the value of s is observed: 
the curves corresponding to different values of s can not be distinguished within 
statistical fluctuations and the fluctuation-dissipation ratio X(t,s) displays a nice 
linear behaviour over a large range of times as can be seen on Figure 1 (on the 
right). To lighten the figures, averages over twenty points have been plotted. On 
the other hand, all data have been used in the fitting procedure. The asymptotic 
value Xoo of the fluctuation-dissipation ratio is obtained by a linear fit of X(t, s) 
with respect to s/t in the range s/t £ [0;0.15]. Again, the data being correlated, 
standard error on the coefficients of the fit underestimates the true error. Taking 
into account autocorrelations of the estimates of X(t, s) leads to an important cor- 
rection to the standard deviation on coefficients given by the linear fit : in the case 
of the Ising model, the linear fit gives an error on of 4.6. 10~ 4 while taking into 
account correlations, it is estimated to be 1.5.10~ 3 . The interpolated line is the bold 
line on Figure 1 and the error bar of our final estimate of Xoo has been put along 
the y-axis. Our final estimates of X^ are collected in Table 1. In the case of the 
4-state Potts model, the fluctuation-dissipation ratio X(t, s) is very noisy. However, 
a rather precise estimation of Xoo is obtained because not all points have large error 
bars so that the main contribution to the fit is due to points with smaller error bars, 
in practice points corresponding to small values of s. Again, the procedure may be 
problematic if corrections depending on s are important. This seems not to be the 
case at regard of the good collapse of curves corresponding to different values of s. 
More problematic is the fact that Figure 1 seems to indicate a downward curvature 
at small values of s/t. The linear interpolation lies outside the error bars for the 
points with the smallest values s/t. Reducing the range of the fit leads indeed to 
smaller and smaller values of X x : 0.459(8) for s £ [0; 0.15], 0.455(8) for s £ [0; 0.1] 
and 0.36(13) for s £ [0;0.05]. Note that equilibrium quantities display logarith- 
mic corrections that may also be present for out-of-equilibrium ones. On the other 
hand, the data being correlated, this curvature may not be a general trend but a 
fluctuation. Since all this is highly hypothetic, we will adopt the safer attitude and 
consider the linear fit inappropriate in this case. 
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X/z X 



Models 



Ising model 0.738(21) 0.328(1) 

3- state clock model 0.844(19) 0.406(1) 

4- state Potts model 0.99(12) 0.459(8)* 



Table 1: Exponent X/z and fluctuation-dissipation ratio for the Ising model, 
the 3-state clock model and the 4-state Potts model, t The value of Xoo given 
for the 4-state Potts model has to be considered carefully since a deviation from a 
purely linear behaviour is observed. 

4 Universality for different lattices 

The procedure detailed in the previous section for the Ising model, the 3-state clock 
model and the 4-state Potts model on square lattices was extended to triangular 
and honeycomb lattices. Note that we only considered regular periodic lattices since 
aperiodic or random lattices may change the universality class at equilibrium. The 
results are presented on Figure 2 for the Ising model, Figure 3 for the 3-state clock 
model and Figure 4 for the 4-state Potts model. The effective exponents X/z and the 
fluctuation-dissipation ratios Xoo are collected in Table 2. Exponents X/z appear 
to be lattice-independent. This statement is in agreement with measurements of 
the dynamical exponent z for the Ising model at T — T c on square, triangular and 
honeycomb lattices [6]. Apart that of the Ising model on triangular lattice, our 
estimates of the fluctuation-dissipation ratios for different lattices are compatible 
within error bars. 

The data for the 4-state Potts model are unfortunately too noisy for the esti- 
mates of X/z to be really useful. The downward curvature of X(t, s) at small values 
of s/t observed on square lattice for this model is also present on triangular lattice 
but not on honeycomb lattice where on the other hand a small upward curvature is 
observed. Despite of these curvatures, all extrapolated values of Xoo are compatible 
within error bars. 

X/z X x 
Models Square Triang. Honeyc. Square Triang. Honeyc. 



Ising 

3- state clock 

4- state Potts 



0.738(21) 
0.844(18) 
0.99(12) 



0.739(22) 
0.845(20) 
0.99(17) 



0.731(17) 
0.844(16) 
0.97(8) 



0.328(1) 
0.406(1) 
0.459(S) t 



0.323(1) 
0.402(3) 
0.460(4)* 



0.328(1) 
0.404(1) 
0.467(21)* 



Table 2: Exponent X/z and fluctuation-dissipation ratio for the Ising model, 
the 3-state clock model and the 4-state Potts model on square, triangular and 
honeycomb lattices, * The value of X^ given for the 4-state Potts model have to be 
considered carefully since in both cases a deviation from a purely linear behaviour 
is observed. 
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Figure 2: Top: effective critical exponent X/z calculated by power-law interpo- 
lation of C(t,s) over the range t € [i m in; 1500] versus s/t m i n for the Ising model 
on the square, triangular and honeycomb lattices (from bottom to top). Bottom: 
fluctuation-dissipation ratio X(t, s) versus s/t for the same models (from left to 
right). The different symbols correspond to s = 10 (circle), s = 20 (square), s = 40 
(diamond), s = 80 (triangle up) and s = 160 (triangle left). The lattice size is 
L = 256. 



5 Universality for different models 

The procedure was also applied to different models belonging to the same univer- 
sality class at equilibrium. We first studied the Ashkin- Teller model at the point of 
its exactly-known critical line where the model belongs to the 3-state Potts model 
universality class. The effective exponent X/z and the fluctuation-dissipation ratio 
Xoo that we numerically obtained are plotted on Figure 5 and our final estimates 
are collected in Table 3. Our estimates of the exponent X/z are incompatible for 
the 3-state clock model and the Ashkin- Teller model. Both are compatible within 
error bars with the estimates found in the literature : 0.828(2) for the 3-state clock 
model [8] and 0.798 for the Ashkin- Teller when fitting the data given by [12] along 
the critical line and assuming a dynamical exponent equal to that of the 3-state 
Potts model [8]. While these models belong to the same universality class at equi- 
librium, it seems thus not to be anymore the case out-of-equilibrium. Surprisingly, 
the fluctuation-dissipation ratios are compatible within error bars for the two 
models, although a small downward curvature may be observed for the Ashkin- Teller 
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Figure 3: Top: effective criticai exponent X/z caicuiated by power-law interpolation 
of C(t,s) over the range t G [t m i n ; 1500] versus s/i m i n for the 3-state clock model 
on the square, triangular and honeycomb lattices (from left to right). Bottom: 
fluctuation-dissipation ratio X(t, s) versus s/t for the same models (from left to 
right). The different symbols correspond to s = 10 (circle), s = 20 (square), s = 40 
(diamond), s = 80 (triangle up) and s = 160 (triangle left). The lattice size is 
L = 256. 



model. 



Models 



X/z 



3- state clock model 0.844(18) 0.406(1) 
Ashkin-Tcllcr model 0.802(20) 0.403(8) 

4- state Potts model 0.99(12) 0.459(8) t 
Baxter-Wu model 1.13(6) 0.548(15)t 
Multispin Ising model 0.977(25) 0.466(3) 



Table 3: Exponent X/z and fluctuation-dissipation ratio Xoo for different models 
belonging either to the 3-state Potts model at equilibrium or the 4-state Potts 
model, t The value of Xoo given for the 4-state Potts model and the Baxter-Wu 
model have to be considered carefully since in both cases a deviation from a purely 
linear behaviour is observed. 
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Figure 4: Top: effective critical exponent X/z calculated by power-law interpolation 
of C(t,s) over the range t G [t m in; 1500] versus s/t m [ n for the 4-state Potts model 
on the square, triangular and honeycomb lattices (from left to right). Bottom: 
fluctuation-dissipation ratio X(t,s) versus s/t for the same models (from left to 
right). The different symbols correspond to s = 10 (circle), s = 20 (square), s = 40 
(diamond), s = 80 (triangle up) and s = 160 (triangle left). The lattice size is 
L = 256. 



We then turned to the study of a multispin Ising model and of the Baxter- Wu 
model both belonging to the 4-state Potts model universality class at equilibrium. 
The effective exponent X/z and the fluctuation-dissipation ratio X^ that we nu- 
merically obtained are plotted on Figure 6 and our final estimates are collected in 
Table 3. Our estimates of X/z for the 4-state Potts model and the Baxter- Wu model 
are compatible within error bars but the latter are very large. On the other hand, 
estimates of X/z for the multispin Ising model and the the Baxter- Wu model are 
fully incompatible. Our values are compatible, although at the boundary of error 
bars, with the estimates found in the literature for the 4-state Potts model [13] : 
0.919 assuming z = 2.294 and the Baxtcr-Wu model [13] : 1.058(4) but not for 
the multispin Ising model [14] : 0.902(10) assuming the same value for z (0.870(11) 
with z — 2.380(4)). Note that these values were calculated from the estimates of 
9 and z obtained by short-time dynamics Monte Carlo simulations. The value of 
X/z that we give is very sensitive to the accuracy of the dynamical exponent z. 
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Figure 5: Top: effective criticai exponent X/z caicuiated by power-law interpolation 
of C(t,s) over the range t G [t m i n ; 1500] versus s/i m i n for the 3-state clock model 
and Ashkin- Teller at the self-dual point y = 3/4 (from bottom to top). Bottom: 
fluctuation-dissipation ratio X(t, s) versus s/t for the same models (from left to 
right). The different symbols correspond to s = 10 (circle), s = 20 (square), s = 40 
(diamond), s = 80 (triangle up) and s = 160 (triangle left). The lattice size is 
L = 256. 



Note as well that we would have obtained smaller estimates of X/z if we would have 
done shorter simulations or restricted our calculations to smaller values of t/s. The 
fluctuation-dissipation ratio reproduces the same tendency as can be seen on 
Table 3. The estimates for the 4-state Potts model and the multispin Ising model 
are in agreement while that of the Baxter- Wu model is clearly not. However, the 
data for both the 4-state Potts model and the Baxter- Wu model show a downward 
curvature as can be seen on Figure 6. The final estimates of have to be taken 
carefully for these models. 



6 Conclusions 

We have addressed the question of universality for ageing ferromagnets, focusing on 
the study of the autocorrelation decay exponent X/z and the fluctuation-dissipation 
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Figure 6: Top: effective criticai exponent X/z caicuiated by power-law interpolation 
of C(t, s) over the range t € [tmin! 1500] versus s/t m ; n for the 4-state Potts model, 
the Baxter- Wu model and a multispin Ising model (from bottom to top). Bottom: 
fluctuation-dissipation ratio X(t, s) versus s/t for the same models (from left to 
right). The different symbols correspond to s = 10 (circle), s = 20 (square), s = 40 
(diamond), s = 80 (triangle up) and s = 160 (triangle left). The lattice size is 
L = 256 for the 4-state Potts model and L = 258 for the others. 



ratio Xrx,. It turns out that, apart from for the Ising model on triangular lat- 
tice, both X/z and do not depend on the particular lattice on which the models 
lives. This supports the idea of an extension of equilibrium universality classes 
to out-of-cquilibrium processes. On the other hand, the Baxter- Wu model and a 
multispin Ising model, both belonging to the 4-state Potts model universality class 
at equilibrium, do not share the same exponent X/z. Moreover, the 3-state Potts 
model and the Ashkin- Teller model, both belonging to the same universality class 
at equilibrium, have different exponents X/z but the same fluctuation-dissipation 
ratio Xoo . This raises the question of the relevant quantities sufficient to charac- 
terise the universality class. Let us recall that we have not studied the influence of 
the dynamics itself on the universal quantities. It would be interesting for example 
to check whether the transition rates Eq. 17 and 18 lead to the same exponents 
and fluctuation-dissipation ratios or not. 
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